Optimal cloning of qubits given by arbitrary axisymmetric distribution on Bloch 

sphere 
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We find an optimal quantum cloning machine, which clones qubits of arbitrary symmetrical dis- 
tribution around the Bloch vector with the highest fidelity. The process is referred to as phase- 
independent cloning in contrast to the standard phase-covariant cloning for which an input qubit 
state is a priori better known. We assume that the information about the input state is encoded in 
an arbitrary axisymmetric distribution (phase function) on the Bloch sphere of the cloned qubits. 
We find analytical expressions describing the optimal cloning transformation and fidelity of the 
clones. As an illustration, we analyze cloning of qubit state described by the von Mises-Fisher and 
Brosseau distributions. Moreover, we show that the optimal phase- independent cloning machine 
can be implemented by modifying the mirror phase-covariant cloning machine for which quantum 
circuits are known. 

PACS numbers: 03.67.-a, 05.30.-d, 42.50.Dv, 



I. INTRODUCTION 

The no-cloning theorem [1| states that no quantum- 
mechanical evolution exists which would transform an 
unknown quantum state jV') according to \^) — >■ |V')|V')- 
This is provided by the linearity of quantum mechan- 
ics. The no-cloning theorem guaranties, e.g., the security 
(privacy) of quantum-communication protocols including 
quantum key distribution. 

Exact cloning is impossible; however, imperfect (opti- 
mal) cloning is possible as it was first shown by Buzek and 
Hillery 2] by designing a 1 — >■ 2 cloning machine, referred 
to as the universal doner (UC). The cloning machine pre- 
pares two approximate copies of an unknown pure qubit 
state. The UC generates two qubit states with the same 
fidelity F = 5/6. Fidelities of the clones to the initial 
pure state are equal Fi = F2. Therefore, the UC is a 
state-independent symmetric doner. 

It was later shown that for the 1 — ?► M UC, a relation 
between the optimum fidelity F of each copy and the 
number M of copies is given by F = (2M + 1)/(3M) [1]. 
Setting M — )■ 00 corresponds to a classical copying pro- 
cess with F = 2/3, which is the best fidelity that one can 
achieve with only classical operations. 

Further works have extended the concept to include 
cloning of qudits, cloning of continuous-variable sys- 
tems or state-dependent cloning (non- universal cloning), 
which can produce clones of a specific set of qubits with 
much higher fidelity than the UC (Jl-llq (see also re- 
views [13, [3 and references therein) . 

The study of the state-dependent cloning machines is 
important because it is often the case that we have some 
a priori information on a given quantum state that we 
want to clone, but we do not know it exactly. Then, 
by employing the available a priori information, we can 
design a cloning machine which outperforms the UC for 
some specific set of qubits. For example, if it is known 
that the qubit is chosen from the equator of the Bloch 
sphere then the so-called phase-covariant doners (PCCs) 



have been designed [3, [ll| , and it was shown to be opti- 
mal providing a higher fidelity than the UC. 

Phase-covariant cloning of qubits has been further ex- 
plored. For example, Fiurasek Q studied the PCCs with 
known expectation value of Pauli's cr^ operator and pro- 
vided two optimal symmetric doners: one for the states 
in the lower and the other for those in the upper hemi- 
sphere of the Bloch sphere. Hu et al. [1J| studied phase- 
independent cloning of qubits uniformly distributed on 
a belt of the Bloch sphere. Bartkiewicz et al. [12| pro- 
vided an optimal cloning transformation, referred to as 
the mirror phase-covariant cloning (MPCC), for qubits 
of known modulus of expectation value of Pauli's cFz op- 
erator. 

Optimal cloning plays a crucial role in, e.g., quantum 
cryptography. Security analyses of the quantum key dis- 
tribution protocols against coherent and incoherent at- 
tacks using quantum doners can be found in Refs. |19l - 
[23j . Optimal phase- independent doners seem to play a 
special role there. One example of the optimal phase- 
independent cloning machines is the PCC, which can be 
used in an optimal attack on the BB84 quantum key dis- 
tribution protocol |19l - [2l| . Another example is the UC, 
which enables an optimal incoherent attack on the six- 
state protocol 0,[2l[2ll. 

Our paper is devoted to phase-independent cloning, 
which refers to cloning of qubits assuming that their 
distribution (symmetrical around the Bloch vector) is 
a priori known. It is related to the optimal state 
estimation [2J] and to phase-independent telecloning 
and tdemapping [16|. Phase- independent cloning in- 
cludes the majority of all known optimal cloning ma- 
chines as special cases. One of the exceptions is the 
phase-dependent cloning machine recently described in 
Ref. [13. 

In this paper, we find that phase-independent cloning 
can be implemented analogously to the MPCC, e.g., in 
linear-optical systems [2g or quantum dots [iJI (see also 
Ref. (23). 



Phase-independent cloning is an example of the opti- 
mal cloning problem being invariant with respect to the 
discrete Weyl-Heisenberg group (see, e.g., Ref. p^). 

We also show here that the phase-independent cloning 
exhibits sudden change in average single-copy fidelity, 
when the symmetry of the system is reduced from U{2) 
to f/(l). 

Optimal cloning also limits the capacity of quantum 
channels. An example is the Pauli channel and Pauli 
cloning machines analyzed by Cerf [28| . Our results could 
be used in a similar analysis for channels that undergo 
phase-independent damping. 

The paper is organized as follows. In Sec. II, we 
present a general transformation describing an optimal 
symmetric 1 — >■ 2 cloning of a qubit. In Sec. Ill, we de- 
scribe the optimal symmetric 1 — > 2 cloning of a qubit 
knowing a priori its axisymmetric distribution on the 
Bloch sphere. This cloning is referred to as the opti- 
mal phase-independent cloning. In Sec. IV, we analyze 
two examples of such cloning of qubits described by the 
von Mises-Fisher and Brosseau distributions. In Sec. V, 
we present, probably the most important result of the 
paper, the optimality proof for the phase-independent 
cloning transformation. In Sec. VI, we describe a quan- 
tum circuit implementing the optimal phase- independent 
cloning. We conclude in Sec. VII. 



II. OPTIMAL SYMMETRIC CLONING OF 
QUBITS 

Suppose we want to clone a set of qubits, for which 
some characteristic point is in the following a priori 
known pure state: 



|V)=cos-|0) 



s*"^ sin -11) 
2' ' 



(1) 



which is parametrized by polar d and azimuthal ip angles 
on the Bloch sphere. We can express all other qubit 
states as 



mo,^)) 



(2) 



where l-tp) is a state orthogonal to \tp), i.e., {tp\'tp} — 0. 
Note that cos 9 is equal to the scalar product of the Bloch 
vectors f — [{ax), (cry), (ctz)] and describes an elevation 
of an arbitrary qubit with respect to the reference qubit 
\ip), whereas </> angle is an azimuth. Please note, however, 
that 9 is measured from lip). 

The quality of the cloning can be described by the 
fidelity of a single clone defined as 



F,(0,0) = (V(0,0)|ft|V'(e,0)), 



(3) 



where pi = Tr^gi (pout) is the reduced density matrix 
of the ith clone (i = 1,2) and © denotes summation 
modulo 2. To obtain an optimal doner for an arbitrary 
distribution of qubits, we use average single-copy fidelity 



F in a way similar to that used for the MPCC [20] . Thus, 
we express F as 

F=\J J 9{0,^)[Fi{9,^)+F2{9,^)]dcos9d(t), 

(4) 
where g(9, (j)) is an arbitrary distribution satisfying the 
normalization condition 

/ / g{9,(l>)dcos9d(l) = l. (5) 



In the special case of 



9[ 



[S{i} - 9) + S{d - 



0)], (6) 



47r sin i3 
our generalized doner reduces to the MPCC. 

The optimal cloning transformation can be found by 
maximizing single-copy fidelity Fi, where « = 1 or 2, av- 
eraged over the distribution g{9, ip) representing the fre- 
quency of occurrence of the considered qubits [12!| . 

Note that distributions g describe our classical knowl- 
edge about qubits to be cloned. Thus, the distributions 
g are classical although cloning is quantum. Probably all 
works on optimal quantum cloning are concerned with 
cloning of pure quantum states. Indeed, a generic quan- 
tum cloning task for pure states can be formulated such 
that the state to be cloned is drawn from a certain ensem- 
ble of pure states with some a priori probability and sent 
to the cloning machine which should produce replicas of 
the input state. The quality of the doner is then usually 
characterized by mean clone fidelity averaged over the 
whole ensemble. The probability distribution character- 
izing the ensemble of states which one wants to clone is 
given by g in our case. So, g is a classical probability 
distribution on a unit sphere. 

An optimal symmetric 1 —> 2 cloning transformation 
must be symmetric with respect to the swap operation on 
clones. So, if we neglect the phase relations, the optimal 
transformation for an arbitrary distribution of qubits can 
be expressed as 

|V')in|00)i,2 -^ COSa+(cOS/3+|V', V')l,2 -Ksin^+IV^, '0)i_2) 

(g)|V')3 + sina+|*+)i,2|V')3, 
|V')in|00)i,2 -^ cosa_(cos/3_|V', V')i,2 -t- sin /3_|-0, -0)1^2) 
®|V')3+sina_|*+)i,2|V')3, (7) 

where |*±) = {\ip,-ip) ± \ip,-ip))/y/2, (-01-0) = and sub- 
scripts 1 and 2 denote the clones whereas 3 denotes the 
ancillary qubit. The angles a± and /3± are found by 
maximizing the average single-copy fidelity. 

It is worth noting that we are not analyzing an optimal 
cloning of mixed qubit states given by 



P 
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g(0,0)|^(0,0))(V(e,0)|dcosM0. (8) 



By contrast, we are analyzing cloning of each of the pure 
qubit states from an ensemble separately. While the dis- 
tribution g{9, (j)) is a weight function for a single-copy 
fidelity of a qubit state \ip{9, 4>)). 



III. OPTIMAL CLONING FOR 
AXISYMMETRIC DISTRIBUTIONS 

In the following, for simplicity, we consider only ax- 
isymmetric distributions. The symmetry of the distri- 
butions means that the cloning transformation is phase- 
independent (does not depend on azimuthal angle 0), 
i.e., distributions of qubits to be cloned are symmetrical 
along an arbitrary jf/;) axis. Moreover, we first postulate 
(as in Ref. [Ij]) that /3± = and then (in Sec. V) we 
prove that such transformation is optimal for cloning. 

When d ^ <j) ~ 0, see Eq. ([T|), all the formulas can 
be rewritten in a computational basis by simply substi- 
tuting \tp) and 1-0) with |0) and |1), respectively. Please 
note that one can proceed also in the opposite direction; 
i.e., one can start assuming that the distributions are 
centered around poles of the Bloch sphere (a local coor- 
dinate system) and later transformed to the basis so that 
the states match in, e.g., the laboratory frame (a global 
coordinate system). We choose the first approach, since 
we believe that it gives better physical intuition of the 
problem. 

The density matrix of a single clone pi is derived by 
taking a partial trace (over the ancillary qubit and one of 
the clones) of the density matrix describing the system 
after the cloning transformation. As a result we obtain: 



in Ref. hi 



6V^ai(a2-l) 



X+X- 



(13) 



where x± — 1 + 2a2 ± 3ai. We notice that for a normal- 
ized qubit distribution (oq = 1) one needs to know only 
ai and 02 in order to fully characterize the corresponding 
optimal cloning transformation. As long as |r| < 1, we 
can express the parameters describing the cloning trans- 
formation as 



2a± — arcsinJl ± arcsinF, 



where 



n^ 



2\/2(l + 2a2)(l-a2) 
■\/3a;+a;_ (3 + ia"^ — iaf — 4a2 



(14) 



(15) 



However, if |r| > 1 then a+ — and a_ = -I or vice 
versa, and the optimal cloning transformations corre- 
sponds to one of the transformation derived by Fiurasek 
in Ref. [9]. The case of Oi = includes the PCC for 
9 = tt/2 and the MPCC [l|] (a+ = a.). Moreover, for 
ai = a2 — 0, we recover the UC transformation [2| . Also 
the optimal cloning transformation of Hu et al. [l3| can 
be derived for an arbitrary belt of the Bloch sphere. 



1 



(cos^ a+ + 1) cos^ - + sin^ a_ sin^ - ) \tp){ip\ 

A A , 



+ ( (cos^ a_ + 1) sin^ - + cos^ a+ cos^ - ) |i/')(V'| 
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sin 6* sin 51 1-0) (-01 + h.c. 



(9) 



where i = 1, 2 and fl = a+ + q;_. Now, we can express 
the single-copy fidelity, given by Eq. ([3]), as 



Fi = - 2(3 + cos2a+)cos'*- + 2 (3 + cos2a_) sin** - 
8 2 2 

-I- [sin^ a+ -I- sin^ a_ + 2\/2 sin ?t\ sin^ Q . (10) 

We obtain expressions for a± by maximizing, with re- 
spect to g(0)-, the average single-copy fidelity F, given 
by Eq. ^^. Any distribution axisymmetric function g(&) 
can be expressed in terms of the Legendre polynomials 
P„ (cos 9) in the following way [23] : 



-. 00 
5W = ^ ^(2n-f l)a„P„(cos0), 



4-77 



(11) 



n=0 
27r /.I 



.9(0)P„(cos6')dcos6'd0. (12) 



Thus, we present the results in terms of the expansion 
coefficients a^. 

The form of the expressions for aj- depends on a single 
parameter P, which is a generalization of the T parameter 



IV. EXAMPLES OF PHASE-INDEPENDENT 
CLONING 



Optimal cloning of the von Mises-Fisher 
distribution 



As an example of phase-independent cloning, let us 
analyze the cloning of qubits described by the von 
Mises-Fishcr distribution (also called the Fisher distri- 
bution) [30] , which is an analog of the normal distribu- 
tion but on a two-dimensional sphere in R (see Fig. [ij. 
The distribution describes dispersion on a sphere and has 
applications not only in physics and mathematical statis- 
tics (especially directional statistics) but also in quanti- 
tative biology, geology, or text data mining. It is a single- 
parameter distribution which could be used to describe 
qubits which undergo a random damping process; for ex- 
ample, weak narrow-band light pulses propagating via 
real media or spin-1/2 systems in the presence of mag- 
netic field fluctuations. The distribution converges to a 
two-dimensional Gaussian distribution for large values of 
K. For K = it becomes a uniform distribution on a 
sphere. 

The von Mises-Fisher distribution is given by 



m 



Kexp (kcos( 
47r sinh k 



(16) 



where k is the concentration parameter (inverted vari- 
ance), which can also be negative when describing dis- 
tribution concentrated around ]0). For the von Mises- 
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(b) 
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FIG. 1. (Color online) Qubit states described on the Bloch 
sphere by the von Mises-Fisher distribution (where the sym- 
metry axis is rotated from the north pole by i9 = 7r/3 and 
</3 — 7r/35) divided by its largest value for (a) k = 0, (b) 
K — 0.3305, (c) K = 1 and (d) k — 10. The darker the region, 
the higher the values of the probability distribution function. 
The mean value (az) = 1. The optimal cloning machines are 
the UC [2] for case (a) and the PCC [y] for the other cases. 
Note that the transition area between (a) and (b) corresponds 
neither to the UC nor the PCC. 
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FIG. 2. (Color online) Average single-copy fidelity F vs con- 
centration K. The dashed line for F — 5/6 corresponds to 
the UC limit, which is reached for k = 0. The shaded area 
corresponds to the range —0.3305 < k < 0.3305, where the 
PCC is no longer optimal. 



on a sphere is the Henyey-Greenstein function [32| . for 
which a„ = /i" and h is the anisotropy factor equal to 
the average of cos9. The Henyey-Greenstein function 
is widely used in simulations of, e.g., volume scattering 
processes. The explicit form of this function is similar 
to the Brosseau distribution discussed in the following 
subsection. 



B. Optimal cloning of the Brosseau distribution 



Fisher distribution, we get the expansion coefBcients, de- 
fined in Eq. p^ . as follows: 



ai = — 1/k + cotliK, 
a2 = 1 — Skai. 



(17) 



If K > 0.3305 then the PCC transformation of 
Fiurasek's Q is optimal. On the other hand, if k = 
the optimal cloning transformation corresponds to the 
UC [3]- So, practically, a very small axisymmetric con- 
centration of qubits (small value of k) , as shown in Fig. [T] 
(Fig. [2]), is enough to determine which of the two cloning 
transformations is optimal. 

We can see in Fig.[5]that the optimal cloning is related 
to the symmetry of the cloned set of qubits (see, e.g., Ref. 



31|). As one would expect, there is a sudden change in 
the average single-copy fidelity when the set of cloned 
qubits reduces its symmetry from U{2) to U{1). 

It is somewhat surprising that the von Mises-Fisher 
distribution, which is an analog of Gaussian distribu- 
tion on a sphere, is not widely used in Monte Carlo 
simulations. This is because of its computational over- 
head. A more popular one-parameter phase distribution 



Another example of phase-independent cloning is 
cloning of photons of known statistics of a single Stokes 
parameter, e.g.. Si. Statistics of the Stokes parameters 
was studied in detail by Barakat .33.]; however, for us 
the statistics of the normalized Stokes parameters is of 
greater interest since it describes the state of a single 
photon. 

In the quantum regime, the normalized Stokes param- 
eters Si — Si/ So for i — 1,2,3 correspond to Pauh's 
operators in the following way: si — Uz, S2 = Ox^ and 
S3 = cTy. Here, the matrices are defined in the basis of 
IIV') = l-ff), \4>) = 1^)}, i-e., for the horizontal and verti- 
cal polarization states, correspondingly. 

For a Gaussian stochastic plane- wave field, the normal- 
ized Stokes s\ parameter has the followi ng p robability 
distribution function derived by Brosseau [3J|: 



m = 



(1 -P2)(l_^c0S( 



2[(1 - ^C0S6')2 - (1 - C0S2 6I)(P2 _ ^2)]t 



(18) 



where cos Q = s\\s, the Stokes parameter, P is the degree 
of polarization, and /i = {Si)/{So); moreover, P"^ — fj? > 
0. The Brosseau distribution has the same asymptotic 
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FIG. 3. (Color online) Average single-clone fidelity F vs pa- 
rameter /i for y? = P^ (solid) and fi ^ P = I (dashed curve). 
The dotted line for F = 5/6 shows the UC limit, which is 
reached for fi = P = 0, whereas the dashed curve corresponds 
to the fidelity of the PCC. 
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FIG. 4. (Color online) Same as Fig. 1, but for the Brosseau 
distribution assuming (a.) n = 0, P — 0.99, (b) n — 0, P — 
0.65, (c) 11 = 0.8, P = 0.99, and (d) ^ = 0.8, P = 0.8. The 
optimal cloning machines are the MPCC Q^ with A — cos a+ 
for case (a) and the PCC 3 ^'^^ ^^e other cases. 



behavior as the von Mises-Fisher distribution if analyzed 
for two hniiting values of P. For P = 0, we have unpolar- 



ized light and the resulting distribution is uniform, which 
corresponds to the UC. For P — 1, the polarization can 
be described as cos6'|i?) -I- exp(i(/)) sin0|T^), where 6 is 
fixed and (f> is an unknown constant. So, the distribution 
converges to Dirac's 6 function 6{cos 9 — fi), which corre- 
sponds to the PCC 0. Then the distribution is given by 
two physical parameters P and fj, describing our knowl- 
edge about the polarization state of the photon. The 
lowest degree of polarization is given as P^ — fj?, this 
corresponds to a situation, in which phase is random 
(see Fig. [3]). By varying those two parameters, we can 
construct various qubit distributions, as shown, e.g., in 
Fig. m For the Brosseau JSJ] distribution, the Ui coeffi- 
cients are more than those in Eq. ^Tl} as given by: 

l-p2 
oi = — ^ — ih - ^.h), 

a2 = \{l-P^){I2-^^h)-\, (19) 



where 



/« = 



1 (l-f M^ - ^^ - 2x/x + a;2p2)3/: 



■da;. (20) 



It is worth noting that the field discussed by Brosseau 
in Ref. [34j was analyzed using the classical description 
of the polarized light. Namely, he considered a Gaus- 
sian stochastic plane-wave field. In his main formula, 
Brosseau used the following notations: the angular brack- 
ets denote the time average, while P is the classical de- 
gree of polarization. So, we use Eq. P^ . which was ob- 
tained using a classical description of light for a quantized 
field. Thus, one can raise an objection that the use of a 
final classical distribution g{9), obtained using the clas- 
sical degree of polarization and the time average of the 
Stokes parameters, in the case of a quantum description 
of the field is incorrect. 

To show that the application of classical distributions 
g{9) for quantum cloning is indeed correct, let us clarify 
the following point: The optimal quantum phase-cloning 
described in the present manuscript and in the vast ma- 
jority of papers on quantum cloning (see reviews [17|, ll8| 
and references therein) is not cloning of mixed qubit 
states, but it is cloning of a pure qubit state from an 
ensemble. We know a priori that the ensemble is de- 
scribed by some distribution g{9) which represents our 
classical partial knowledge about the qubit to be cloned. 
Thus, the distribution is classical although cloning trans- 
formation is quantum. 

Our application of classical distribution g(9) is in com- 
plete agreement with other works on optimal doners (see 
Table HI . Our examples include (i) the universal doners 
of Buzek and Hillery p] , where the classical distribution 
g{9) is equal to j-, (ii) the optimal phase-covariant don- 
ers of Brufi et al. %[^ (for ^ = 7r/2) and of Fiurasek [10 
(for any i?), where g{9) is given by Dirac's S function 
2Ti-sin(i?') '^(^ ~ ^)' (^^^) ^^^ optimal mirror phase-covariant 
doners of Bartkiewicz et al. [12j, where g{9) is given by 
Eq. (6), and (iv) the optimal quantum doners of a state 



in a belt of Bloch sphere as described by Hu et al. [l3|, 
where g{0) is given by a rectangle distribution, i.e., states 
are assumed to be uniformly distributed in a region form- 
ing a belt between two latitudes on the Bloch sphere. 



V. OPTIMALITY PROOF FOR THE 

PHASE-INDEPENDENT CLONING 

TRANSFORMATION 

Cloning transformations can be described by a com- 
pletely positive trace-preserving map (CPTP) x^ which 
can be written with the use of the inverse Jamiolkowski 
isomorphism [S^ as follows: 



Pout = Tri„(xp4®l^') 



(21) 



The average single-copy fidelity (the figure of merit when 
maximized) can be expressed as 



i^-Tr(xi?), 



(22) 



and 



^=2 



2ir 



pJn ® (Pin «) 1 -t- 1 «) Pin) g{e) d COS Acj). 

(23) 
Since the cloning transformation for the set of cloned 
qubits needs to be symmetric with respect to an arbitrary 
rotation U about the state \'ip) and to swapping C/swap 
of the clones, the following commutation relations should 
hold: 



[x, u* ® U^^] = 0, 
[x, a ® UswAp] = 0. 



(24) 



The commutation relations and trace-preserving condi- 
tion restrict the x rnatrix to the following form (in 

{|V'),|^)}®3 basis): 
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(25) 



where ^4 = 1 — 2^2 — Ci ^-nd 7/4 = 1 — 2?72 — ?/i. It is 
more convenient to use another basis to express x and 
i? operators, i.e.: {IV^,-;/^, V)^|Vi)|*+), IV^,^, V-), IV')I*+), 
|f/')|*-), |-0)|*-),|V',V','0),|V')V')V')}- In this new basis 
the operators have the following block form: 



0», R^^R^. 



(26) 



i=l 



where Xi(2) and Ri{2) are 2x2 matrices. The Xi and Ri 
elements for i = 3, ..., 6 are not negative numbers limited 



by trace of x (Tr x = 2). The average single-copy fidelity 
can be expressed using Einstein's summation convention 
as 



F = Tr(i?,x.)- 



(27) 



The fidelity is a convex superposition and the matrix 
elements Xi (for * = 3,4,5, and 6) are numbers € [0,2]. 
Therefore, it can be proved that for any distribution g{9) 
for i e {5,6} we have Ti{RiXi) < 1/2. Moreover, for 
i e {3,4} we have Tr (i?iXi) < -Fpcc, where Fpcc is 
the single-copy fidelity of the PCC averaged over inputs 
characterized by g{9). The PCC is derived only from xi 
and X2- 

Hence, in order to maximize the fidelity we have to 
put X3 = X4 = X5 = X6 = 0. This implies that r/2 = %, 
^2 = £.3 and 1 - ?7i = 2772, 1 - Ci = 2^2- Now, the 
CPTP map depends only on four parameters and can be 
described as a direct sum (denoted by 0) of two matrices: 



X 



m V2C1 
\/2Ci 1-6 



Ci V2C2 
72(2 i-m 



(28) 



where the first matrix acts in a subspace 
spamied by {\ip , ^p , ip) , \^ , "$ +)} and the second by 
{iV'j "01 "0)7 IV'i ^-i-)}- The map can be parametrized 
without the loss of generality in the following way: 



X 



72(1 



V2C1 

• 2 
sm a_ 



V2C2 



72(2 

• 2 

sm a+ 



(29) 



However, for the extremal CPTP maps [3i[, we have 
that Xi(-2~i ^ ^1(2)- From this condition, we find that 

\/2Ci = cosa+ sina_ and a/2C2 = cosa_ sina_(_. Finally, 
any CPTP map that maximizes the average single-copy 
fidelity for an arbitrary axisymmetric distribution g(9) 
can be written in the basis {10), IV')}'^'^ as 
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(30) 



\/2 V2 

where c± = cos a± and s± = sin a± . The above map 
can be decomposed into unitary transformations, given 
by Eq. Q ( for /3± — 0), by means of the Kraus decom- 
position [36| . This completes the proof. 



VI. QUANTUM CIRCUIT FOR OPTIMAL 
PHASE-INDEPENDENT CLONING 



The analyzed cloning problem can be also expressed 
in the logical basis. The optimal cloning transformation 



TABLE I. Summary of optimal phase independent symmetric 1 — >■ 2 cloning machines for qubits. 



Optimal cloning machine 



Distribution g{8) 

X 

47r 

2tt sin -I? 

5{-&-e)+5{-&+Tv-e) 

An sin i? 

u{e-^i)-u{e-^2) 

27r(cos i?i —cos i?2 ) 

K exp (k cos 6) 

4tt sinh k 

(l-P^)(l-^cos6') 



Single-copy fidelity F 
5 



I [5 + ^/2 + 2 cos(6' + k) - (^ - 1) cos(26i)] 



Universal doner (UC) 

Phase-covariant doner (PCC^ 

Mirror phase-covariant doner (MPCC)^ 

Cloner of Bloch-sphere belt_^'_^ 

Cloner of von Mises-Fisher distribution^^ 

Cloner of Brosseau distribution*^ 



l+A^ 



I (A^ - AAV2) 



2[(l-Atcos6»)2-(l-cos2 e)(p2-^2)]3 



2 2 

Eq. ^ with Eq. (15^^ form Ref. [JJ 
Eq. ^ with Eq. ^ 
Eq. ^ with Eq. ^ 



K = for < i9 < I and K = TT for I < i9 < TT. 



''A=y1— A2,A depends on i? as given in Ref. [12| | . 

•^ The explicit expression for F is lengthy and therefore is not shown here, but it can be directly derived form the quoted equations. 

u is Heviside's step function. 
" In Ref. |13|| a = 0+ and [3 = o_. 



can be now written as 

|000) ^-cosa+lOOl) -f sina+|V^+)|0), 

|100> ^cosa_|110> +sina_|^+)|l), (31) 

where |V'+) = (|01) -(- |10)) /v^. The quantum circuit, 
shown in Fig. [SJ performs the fohowing transformation: 

liPout) — I^CNOT'^CNOT'^CNOT'^CH "^Ry '^CRyl'f^in/' 

(32) 
where the superscripts indicate qubits for which the cor- 
responding gate is applied. The basic elements of the 
circuit are the rotation 



C/r,M- 



cos (w/2) — sin (w/2) 
sin (a;/2) cos (a;/2) 



(33) 



about the y axis by angle a; — loij^ and the controlled 
rotation ?7cRy('&) by angle $ — 2(a_ — q+). In ad- 
dition, this circuit is composed of the controlled NOT 
(CNOT) gates, C/cnot, and the controlled Hadamard 
gate, C/cH, which can be decomposed [12j as C^ch — 



^ "-^CNOT^ ' 



where 



\f^^^^2^ 



1 l + \/2 
l + \/2 -1 



(34) 



|0). 



R 



|0)3-i?y($) 



i?j,(a;) 



^^ e 



■^ 



^— ^ 



FIG. 5. A quantum circuit which implements the opti- 
mal phase-independent cloning transformation of |i/;(S,0)) = 
cos (61/2) 1 0) + exp (z(^) sin (61/2) 1 1) described on the Bloch 
sphere by arbitrary axisymmetric distribution function g(Q\ 



From left to right: 
by angle $ = 2(a_ 



controlled rotation Ry about the y axis 
— a+), rotation Ry by angle oj = 2a+, 



controlled Hadamard gate, and three CNOT gates. 




FIG. 6. A quantum circuit, shown in Fig.O implementing the 
optimal phase-independent cloning transformation, is equiv- 
alent to the MPCC [12] with A = cos a+ together with the 
controlled rotation Ry by angle $ — 2(q_ — a+)- 



VII. CONCLUSIONS 



The optimal phase-independent cloning can be imple- 
mented with the use of a quantum circuit, the MPCC [12] 
setting for A — cosa+ (see Fig. [S]). Therefore, with mi- 
nor modifications (concerning the state preparation of 
the third qubit), the optimal phase- independent cloning 
machine can be realized with, e.g., linear optics orquan- 
tum dots as described for the MPCC in Refs. [ll,[2l[33l. 



We analyzed optimal state-dependent cloning of qubit 
states, which are described by a priori known arbitrary 
phase-independent (axisymmetric) distribution g on the 
Bloch sphere. This optimal cloning reduces in special 
cases to the universal cloning of Buzek and Hillery [3], 
the phase-covariant cloning of BruJ3 et al. [5| and its gen- 
eralization by Fiurasek [8, 9] , the mirror phase covariant 



cloning (MPCC) of Bartkiewicz et al. [T2], or cloning of 
an uniform belt of the Bloch sphere of Huei al. |13l |. 

As an example of the state-dependent cloning, we stud- 
ied the cloning transformations of qubits described on 
the Bloch sphere by the von Mises-Fisher and Brosseau 
distributions, where the first is an analog of normal dis- 
tribution on a sphere [30| and the latter describes statis- 
tics of the Stokes parameters [33i, 34] . Whereas the first 
example is more abstract and describes Gaussian-like dis- 
persion, the second example can be used directly to es- 
timate the upper bound for the capacity of a depolar- 
izing channel [28j for photons. Our results can be also 
applied in security analysis of various quantum commu- 
nication protocols, including quantum teleportation [38| 



and quantum key distribution [19J. 

Recently, it was shown that phase independent-cloning 
can be parameterized by four parameters [23 ■ Here, we 
proved that only two parameters are sufficient to describe 
the optimal phase-independent cloning. 

Moreover, we showed that the phase-independent 
cloning is a simple generalization of the MPCC and, 
thus, can be implemented analogously to the MPCC 
using photon-polarization qubits i26lj and quantum-dot 
spins 12]. 
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